UWThPh-2000-26 



A neutrino mass matrix with seesaw mechanism 
and two-loop mass sphtting 

W. Grimus 

Universitat Wien, Institut fiir Theoretische Physik 
Boltzmanngasse 5, A-1090 Wien, Austria 

L. Lavoura 

Universidade Tecnica de Lisboa 
Centre de Fisica das Interacgoes Fundamentais 
Institute Superior Tecnico, 1049-001 Lisboa, Portugal 

3 July 2000 

Abstract 

We propose a model which uses the seesaw mechanism and the lepton number 
L = Lf. — — Lt to achieve the neutrino mass spectrum mi = m2 and = 0, 
together with a lepton mixing matrix U with Uez = 0. In this way, we accommodate 
atmospheric neutrino oscillations. A small mass splitting mi > m2 is generated by 
breaking L spontaneously and using Babu's two-loop mechanism. This allows us 
to incorporate "just so" solar- neutrino oscillations with maximal mixing into the 
model. The resulting mass matrix has three parameters only, since L breaking 
leads exclusively to a non-zero ee matrix element. 

The recent results of Super-Kamiokande [0], providing evidence for atmospheric neu- 
trino oscillations, have lead to increased efforts to investigate mechanisms for the 
generation of neutrino masses and mixings. In this context one wants to solve two ques- 
tions: 

1. Why are the neutrino masses so much smaller than the charged-lepton masses? 

2. How can the specific features of the neutrino mass spectrum and of the lepton 
mixing matrix , needed to reproduce the atmospheric and solar- neutrino deficits, 
be generated within a model? 
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Proposals to answer the first question are given by the seesaw mechanism ||^ and by 
radiative neutrino-mass generation. If we confine ourselves to extensions of the Standard 
Model in the Higgs sector then purely radiative neutrino masses are obtained within 
the models of Zee [|], 0] and of Babu and extended versions thereof |^ 0, |Tl|. As 
for the second question, one prominent feature of the mixing matrix U is the smallness of 
for which one would like to find an explanation in a model. As a means to 



Ue3 112, 13f 



Ln — Lj. has been suggested [14 



achieve this, the lepton number L 

In this letter we propose a model which combines the seesaw mechanism, Babu's 
radiative two-loop mechanism, and the lepton number L. In this way, we have an ex- 
planation for the smallness of the neutrino masses. The seesaw mechanism will enable 
us to fit the atmospheric neutrino oscillations, whereas Babu's mechanism will generate 
the small mass-squared difference necessary for solar-neutrino oscillations. The lepton 
number L will insure f/gs = 0. However, the breaking of this lepton number is crucial for 
solar-neutrino oscillations. 

Our model is given by the Standard Model of electroweak interactions, based on the 

gauge group SU(2)xU(l), with three Higgs doublets (pk = (v^fc^V^fc) ■, where k = 1,2,3. 
The vacuum expectation values of V^fl are denoted Vk- We also introduce two neutrino 
singlets ujij (j = 1, 2) and the scalar singlets h'^^, and rj. The latter scalar is complex 
but has zero electric charge. We have the following assignments of the lepton number L 
to these multiplets: 
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Furthermore, we need a discrete symmetry S defined by 



S : Hr {i = e,fi,T), 



-^02, 



V -V, (2) 



while all other multiplets, in particular the left-handed lepton doublets, transform trivially 
under S. With Eqs. (|I|) and we obtain the Yukawa Lagangian 
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J2 rUiiuiLji) 
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0* 



K i'^eL, eL)h'Rl + 5* (Z/^L, /iL)^'iJ2 + K i'^rL, TL)h'R2 

+ E hu'fRC-H'Rh++ + 



(3) 



where 5a (a = e,/i, r), fr, and the are complex coupling constants. The Yukawa 
Lagrangian in Eq. @ is the most general one built out of the multiplets in our model and 
compatible with the gauge symmetry and with the symmetries L and S. The first line of 
Eq. displays the ordinary Yukawa couplings to the Higgs doublet 02, which give mass to 
the charged leptons. Notice that we have taken, without loss of generahty, those Yukawa 
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couplings to be flavor-diagonal. The second line of Eq. shows the Yukawa couplings 
of the right-handed neutrino singlets to 0i; the third and fourth line display the Yukawa 
couplings needed to implement radiative neutrino masses using Babu's mechanism 
Notice that the symmetry S forbids the couplings u^-^C~^URirj and h']^2C'~^^R2V* j which 
would be allowed by L. 

Since the right-handed neutrino fields are gauge singlets, we can write down the mass 
term 

CM = M*ul,C-'iyn2 + il.c., (4) 

which is compatible with the lepton number L. In the following, M will play the role of 
a large seesaw scale. 

The neutral complex scalar rj breaks the lepton number L spontaneously through its 
vacuum expectation value (?7)o- In this way, the only term in the Higgs potential which 
is linear in r], 

V, = Xvf-f-h+^ + }l.c., (5) 

where A is a dimensionless coupling constant, transforms upon spontaneous symmetry 
breaking into 

^sb = A(r7)orr/i++ + H.c. (6) 

V^b provides the trilinear scalar coupling required by Babu's mechanism IQ. Note that 
one cannot introduce a priori the Z-soft-breaking term f^h~^^, since this term has 
dimension 3 and then we would also have to introduce the mass terms z/^jC~^z//jj (j = 
1, 2), which also break L softly and have dimension 3; else we would not have a technically 
natural model . These mass terms would destroy the L invariance of the light-neutrino 
mass matrix already at tree level. 

The symmetry S necessitates the introduction of two Higgs doublets in order to have 
enough freedom to make the charged leptons massive. The Higgs doublet 03 does not 
couple to the leptons due to the symmetry 5", but it is needed in order to have the terms 
(0103)^ and (0102) (0302) in the Higgs potential. These terms prevent the appearance of 
a Goldstone boson which would couple directly to the leptons. On the other hand, the 
spontaneous breaking of L results in a Majoron, given, if we assume that (77)0 is real, 
by \/2lmr]. However, this Majoron is only very weakly coupled to matter (electrons, up 
quarks, and down quarks) via loop diagrams 



Notice that, due to the specific form of the symmetry S, there are no couplings of 
the type /~0fcT20A;', where T2 is the antisymmetric Pauli matrix and therefore k 7^ k'. 
The absence of these couplings impedes Zee's mechanism for one-loop neutrino masses 
0, contrary to what happens in other models [Q. 

Let us now discuss the neutrino mass matrix. We have a 5 x 5 Majorana mass matrix 
Ai following the five chiral neutrino fields in our model. The neutrino mass term is given 
by 

£d+M = ^ nlC-^MQL + H.C, with Ql = {l^eL, y^,L, l^rL, (z^Rl)', (^^ffi)')^ • (7) 

The mass matrix has the decomposition 
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with, at tree level, 



according to the second line of Eq. and to Eq. (^. The matrix vanishes at tree 
level. Thus, at tree level we have [1^ two degenerate neutrinos with large masses, two 




degenerate neutrinos which are light due to the seesaw mechanism, and one massless Weyl 
neutrino. But, since the lepton number L is broken by V^b in Eq. (H), Ml 7^ is generated 
at two loops by Babu's mechanism. However, as L is operative in the Yukawa couplings 
of and of h^~^, one finds that only (Mi)ee is non-zero. We obtain 

Ml= \ I , with TUee = ^^^^^ J2 fehu>femfmfJu>. (10) 

Here, Iw is a convergent two-loop integral. Assuming that the masses of and of h'^'^ 
are of the same order of magnitude and are much larger than the masses of the charged 
leptons, one finds ||, [1^ /«/ ~ l/m|, since all (di-)logarithms are of order 1. 

We may now apply the seesaw formula in order to obtain the neutrino Majorana mass 
term for the light neutrinos as 

-^m = ^UjIC~'^M^ UJl + H.C., with UJl = {l^eL, T^f,L, l^rhf , (H) 

where 

M^ = Ml- MIM]^^Md. (12) 
Inserting the expressions for Ml, Md, and Mr into Eq. (|13), we arrive at 

/ rrtee ~SJ^/M -6Jr/M \ 

= -5e5jM • (13) 

V -5e5r/M ; 

Obviously, the second term on the right-hand side of Eq. (|r^), i.e., Eq. (|I^) with rriee = 0, 
is invariant under the lepton number L. 

By phase transformations, all elements of the mass matrix in Eq. (|l^) can be made real 
and non-negative. Consequently, there is no CP violation associated with this neutrino 
mass matrix. The matrix is exceedingly simple, since it is parametrized by only three 
positive quantities a, r, and 

Ml= r6 I . (14) 




^As a matter of fact, the mass matrix in Eq. (|l4| ) is the one imphcitly suggested by Joshipura and 
Rindani at the end of their paper. However, these authors reUcd on Zee's mechanism instead of relying 
on the seesaw mechanism, and consequently they ran into difficulties with the orders of magnitude of the 
various terms needed in order to fit both the atmospheric and the solar-neutrino oscillations. They have, 
therefore, discarded it. 
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The eigenvalues of M^, expressed as neutrino masses, are mi, — m2, and = 0, with 

(15) 



mi 



where we have defined 



'mo + — + - and m2 



/ „ a 

nin H 

° 4 2 



mo 



6v/X, with X = 1 + . 



(16) 



Furthermore, from we derive the lepton mixing matrix 



/ 
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with c = cos 00, s = sin 6', 



70, 



(17) 



where Oq is the solar-neutrino mixing angle. Notice that t/gs = 0. This is an important 
and exact prediction of our model. 

Since the matrix element a in M'^ is meant to generate a small mass splitting between 
mi and m2, necessary in order to accommodate solar- neutrino oscillations, this element 
must be tiny. Therefore, using the neutrino masses in Eq. (15) together with = 0, we 
find 

,2 _2 A™2 2moa (1^ 



Am^^jn — "^0 Am| ~ 



for the atmospheric and solar mass-squared differences, respectively. This in turn gives 

and sin^ 26q ~ 



a 



Ami, 



1 



1 



16 \ Am2 



(19) 



atm , 



As a consequence, the solar-neutrino mixing angle is for all practical purposes 45°, and a 
fit to the solar-neutrino data requires the "just so" or vacuum-oscillation (VO) solution 
to the solar- neutrino deficit (for recent analyses see Refs. |T8| , |20| , |2T|). The so-called 
LOW solution |l^, ^ might also be allowed ||23|. In the following we shall however 
concentrate on the VO solution. Due to Ues = 0, solar-neutrino oscillations are decou- 
pled from the atmospheric oscillation parameters . The mixing angle for atmospheric 
neutrino oscillations is obtained from the mixing matrix in Eq. (|T^) as 



sin^ 26, 



4^2 



:i+r 



2^2 ■ 



(20) 



According to the Super-Kamiokande results |2^, sin^ 26'atm ~ 0.88 at 90% CL, and one 
therefore has the range 

0.7<r<1.4 (21) 



for the parameter r. We then have a scenario very close to bimaximal mixing ^6 
Thus, with Eqs. (||), (0), and 



'Am^ 
1 +r2 ' 



(22) 
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all three parameters of are in principle determined by the atmospheric and solar 
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neutrino oscillation data. Summarizing, with Am%. ~ 10 eV^, Am^^j^^ ~ 3.2 x 10 ^ 



eV [25|, and the atmospheric neutrino mixing angle close to 45°, we arrive at the order- 



of-magnitude estimates 

a ~ 10"^ eV, b ~ 0.04 eV, r ~ 1. (23) 

We now want to study how to implement these estimates in our model. To this end 
it is useful to consider the different mass scales present in the model. Apart from the 
electroweak scale, and with the following simplifying assumptions, we have three scales: 

• The scale, denoted mn, of 6e, S^, and 6r — assuming that these parameters are all 
of the same order of magnitude.0 

• The mass scale of the new scalars, where we make the reasonable assumption ~ 
ruf ~ (77)0. 

• The mass M of the right-handed neutrino singlets. 

Obviously, the two mass-squared differences Aml^^ and Am|, cannot determine all three 



mass scales. However, the two-loop expression in Eq. (10) results in 



1 , „o. 



a ~ 



(167r 



2^2 



XfXr — , (24) 
rrih 



where we have assumed that all the couplings hui are of similar magnitude (and also 
l/^l ^ l/^l) and that, therefore, the two-r contribution to Babu's two-loop diagram is 
dominant. Equation (E3) allows us to estimate the mass scale of the Higgs scalars: 



\f%r X 10^^ GeV. (25) 



With the dimensionless coupling constants A, and h^r being at most of order 1, we 
may consider 10^^ GeV as an upper bound for rrih- Clearly, with the scalar singlets having 
such large masses, there are no restrictions on their Yukawa interactions stemming from 
decays and scattering data. 

From the atmospheric mass-squared difference we get 



M 



A^atm ~ 0.06 eV. (26) 



We can tentatively fix m^) by making the assumption that it is of the order of m^, the 
largest of the charged-lepton masses. This leads to M ~ 10^^ GeV. Interestingly, this 
order of magnitude is compatible with mn — see Eq. (^) — if we assume that the 
dimensionless coupling constants are ~ 10^^. An attractive option would be to identify 
the two mass scales, i.e., to assume M ~ 111^. 

To summarize, in this paper we have advocated the three-parameter neutrino mass 
matrix in Eq. (0). That mass matrix yields maximal solar- neutrino mixing and, there- 
fore, it requires the vacuum-oscillation solution to the solar-neutrino deficit (the LOW 



^(5^ and 5r certainly are of similar magnitude, as follows from Eqs. (fsh, dlj), and pT 
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solution might also be an option). On the other hand, in order to obtain nearly maximal 
atmospheric neutrino mixing one has to tune the parameter r to be close to 1. With this 
mass matrix one gets the neutrino mass spectrum mi ~ m2 and = 0, while 11^3 = in 
the lepton mixing matrix of Eq. (]T7|). Notice that ms = 0, the inverted mass hierarchy, 
and Ue3 = are exact, testable predictions. We have furthermore shown that the mass 
matrix in Eq. (0) can be obtained in an extension of the Standard Model with a sponta- 
neously broken lepton-number symmetry L — the ensuing Majoron couples very weakly 
to matter — together with a seesaw mechanism, responsible for the smallness of the mass 
b, and a radiative mass-generation mechanism, responsible for the tiny mass a. We have 
shown that it is sufficient to have a single heavy mass scale in our model, comprising both 
the seesaw scale and the scale associated with the masses of the new gauge-singlet scalar 
particles; a value of order 10^^ GeV would be a suitable choice for this heavy mass scale. 
The smallness of the mass a in the mass matrix of Eq. practically forbids neutrinoless 
PP decay. 
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